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Abstract 

We investigate a family of integrals involving modified Bessel functions that arise in 
the context of neutrino scattering. Recursive formulas are derived for evaluating these 
integrals and their asymptotic expansions are computed. We prove in certain cases 
that the asymptotic expansion yields the exact result after a finite number of terms. In 
each of these cases we derive a formula that bounds the order at which the expansion 
terminates. The method of calculation developed in this paper is applicable to similar 
families of integrals that involve Bessel or modified Bessel functions. 


1 Introduction 

In this paper we study linear combinations of the family of integrals 

J f' oo 

cosh N (t)K I/ (zcosh.(i))dt, %\(z) > 0, (1.1) 

o 

where the modified Bessel function of the second kind for Dap) > 0 and v e C is given by the 
integral 

r oo 

K u (z) - / exp(-zcosh(£)) cosh.(ut)dt. (1.2) 

Jo 

See, for example, [I] for a detailed treatment of the theory of Bessel functions. 

The following special case is of particular interest 

r oo 

Gp(z)= / cosh N (t)K u (zcosh(t))smh 2 (t)dt, Dip) > 0. (1.3) 

Jo 

Note that G^ can be obtained from the by the linear combination 

G»(z) = F?+\z)-F?(z). (1.4) 

Integrals of the form Eq. (II. lj) appear in the computation of the reaction rate for electron- 
positron annihilation into neutrinos, as outlined in Appendix |A] 
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The integral Eq. (II.ip is related to the Meijer G-function function and Mathematica [2] 
is capable of expressing the general case Eq. (II.ip in terms of generalized hypergeometric 
functions. In addition, a related integral formula is well known, see equation 13-72 (1) in [I] 
or equation 2.16.13.2 in [3]. 

Lemma 1.1. Let fH(z) > 0 and a,b e C. Then 

J f' oo 

K a+ b(2zcosh(t)) cosh((a - b)t)dt - K a (z)K b (z). (1.5) 

o 

By expanding 

N 

cosh N (t) = Y^CkCOsh(kt), (1.6) 

k=0 

Eq. (II.5p can be used to express F^(z) as a linear combination of products of modihed Bessel 
functions. 

In this paper an alternative, and in our opinion simpler, method for computing Eq. (II.Hi 
will be given that also provides more detailed information about the structural form of the 
result in some cases. In particular, we show how to compute these integrals recursively, 
using only the case a - b from Eq. (1 1.5 j) and without explicitly performing the expansion 
Eq. (11.61) . We will also prove that if N is even and v is odd or vice versa then F/f(z) can 
be expressed as a finite linear combination of terms of the form e~ z z~ n , and compute the 
order n at which these expansions terminate. In general, the order is smaller than one would 
naively expect from combining the expansion Eq. (II,6j) with Eq. (II,5j) . More specifically, the 
following theorem will be proven. 

Theorem 1.2. Let N, v e {0,1,2...}. If N is even and v is odd or vice versa then 


where 


F„(z) = z 1 e 


Mn,v 

E 

s=0 


E «,r(r + i/2) y . bk,NClj 


j+r=s 


k+l=r 


jis - 1 if N < v, 
jlV-1 if N > v 


(1.7) 


( 1 . 8 ) 


and 


Kn,i 

G?(z) = z~ 2 e~ z £ 


s =0 Lj+t=s 


r(r + 3/2)y dk,NCl K 


k+l=r 


(1.9) 


where I< N ^ = M N+2>v - 1. 

See Eq. H2.191) through Eq. \ 2. 221) for the definition of the coefficients aj, bj tN , dj^, and 

c Ij- 


Section [2] contains necessary background material on asymptotic expansions, including 
explicit formulas for the asymptotic expansions of Fff(z) and G{f (z) as z -> oo. Section 
[3] contains proofs of several recursion relations that the family of integrals Eq. (11.11) obeys. 
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Using these recursion relations, a proof of Theorem 1 1.21 will be given in section [4j The case of 
interest to neutrino scattering is N and v both odd. Although Theorem 11.21 does not apply, 
the recursion relations derived in section [3] still allows for an analysis of this case. This is 
discussed in Appendix IA1 


2 Asymptotics 


In this section we collect some needed results on asymptotic expansions of integrals with ex¬ 
ponentially decaying integrands. For these purposes, it is convenient to rewrite the integrals 
Eq. (11.11) and Eq. (11.31) as 

F v( z ) = f y T + 1 ) K v (z(x + l))dx, (2.1) 

~'° \Jx{x + 2) 

Gf (z) ~ [ \/ 'x(x + 2)(x + 1) N K u (z(x + 1 ))dx. (2.2) 

Jo 

We will consider these integrals for z, N, v e M as z -*■ oo. 

Laplace’s method will be needed. See, for example, pt]. 

Lemma 2.1. Let cj : (0, oo) M and b > 0 such that q e L 1 (e~ zx dx) for all z>b. If 


q ~ E a n x ( ' n+x as x ^ 0 

n=0 


(2.3) 


for A, n > 0 then 


r 

Jo 


l q(x)dx 


£r(^) 


Cl ri 


fl ) 2 (( n+A )/^ 


as 2 : 


00 . 


(2.4) 


In fact, we will also need a generalization of Laplace’s method to integrals where the 
kernel e~ x is replaced by another exponentially decaying function. The general method for 
deriving an asymptotic expansion in this scenario is given by Olver in Section 9.6 of |3]. For 
the present purposes, the following variant is needed. We include a proof for completeness. 


Lemma 2.2. Let K : (0, 00 ) M have an asymptotic expansion 


K(x) ~ x a e ,8x ^ ^ as x 00 

n=0 


(2.5) 


where a e M. (3 > 0 and f e ©(e 0 "*) as x ^ 00 for some a > 0. If K is bounded on compact 
subsets of (0, 00 ) (but not necessarily around 0) and f is integrable on bounded sets then for 
any M > 0, and sufficiently large z 


J f" 00 

f(x)K(z(x +1 ))dx 
0 

r 00 

e ~^ z / f(x)(x + l) a ~ n e~ f8zx dx + 0( 

Jo 


M -1 

n=0 


z a-M e -pzy 


( 2 . 6 ) 


If f is also bounded on bounded sets then the error term can be improved to 0(z a M 1 e I z ). 
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Proof. In this proof we employ the notational device of using the same symbol for different 
constants used in bounding asymptotic behavior, here denoted by C. 

First it must be shown that for sufficiently large z, both sides are defined. For the left 
hand side, fix a > 0 large enough so that |/| < Ce ax and \K\ < Cx a e~@ x on (a, oo) and take 
z>a,a/f3. Then 

J f' oo 

\f(x)K(z(x + l))\dx (2.7) 

o 

J r a r oo 

\f(x)K(z(x + l))\dx + C / e~^ z ~ a>,x ~P z {z{x + l)) a dx < oo. 

0 Ja 

The right hand side is similar. Therefore all integrals are defined for z sufficiently large. 
Now compute 


J r' oo Jv-i r' oo 

f{x)K{z{x + 1 ))dx - Y a n z a - n e~P z / f(x)(x + l) a ~ n e^ zx dx 
o Jo 

/-oo N -1 _ 

< / \f(x)A N (z(x + l))\dx, A N (x) = K(x) - x a e~ /3x E — 

J 0 n-0 % 


( 2 . 8 ) 


Take a sufficiently large so that |A w (x)| < Cx a ~ N e~h x and \f(x)\ < Ce ax on (a, oo) and let 
z>a,cr//3. Then 


r oo 

/ \f(x)A N (z(x + l))\dx (2.9) 

Jo 

<Cz a ^ N e~^ z (J~ \f(x)\(x + l) a ~ N e~ l3zx dx + J (x + l) a ~ N e~^ z ~ a ^ x dx\ 


<Cz a ~ N e~P z . 


If / is also bounded on bounded sets then for any e > 0 

J r oo 

\f(x)A N (z(x + l))\dx 
o 

/ r a r oo 

<Cz a ~ N e~ pz ( / e~ fizx dx + / (x + l) a ~ N e~ (/Sz ~ a)x dx 

\J 0 Ja 

<Cz a ~ N e~ pz (~ + e-^ z ~ a - e)x dx^ 

<Cz a ~ N e~P z ( - + -— - ) 

\z pz - a - e) 

<Cz a ~ N ~ l e~P z 
for z sufficiently large. 

Finally, we need to be able to combine asymptotic expansions. 
Lemma 2.3. Let f,g : (0, oo) -* M such that 


oo oo 

/ ~ E a n x (n+x ~^^, g ~ E b nX t n +h-AA asx-+ 0 

n =0 n-0 


( 2 . 10 ) 


□ 


( 2 . 1 !) 
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for n,\,/3>0. Then 


/s~y E 


( 2 . 12 ) 


r =0 \m+n=r 

Together, the above lemmas give the following theorem. 


Theorem 2.4. Let K : (0,oo) -> M be bounded on compact subsets of (0,oo) and have the 
asymptotic expansion 


I< ~ x a e~' 3x £ 


as x -> oo 


—* rp I i 

-n 


(2.13) 


n=0 


where ael, /3 > 0. Let / e 0(e CT:r ) as x ^ oo for some a > 0, &e integrable on bounded sets, 
and have the asymptotic expansion 


f ~ Y b n x n+X 1 as x -> 0 

n =0 


for A > 0. Then for any M > 0 and z sufficiently large, 

J r oo 

f(x)K(z(x + l))dx 
o 


M-l 

= E 


E 


ajT(r + A) 


s=0 Lj+t=s 


p 


r+A 


Y bkCl,j,a 


k+l-r 


^a-X-s e ~f3z + Q( z a-\-M e -pzj 


where 


(x + l) a 3 ~ Y c i,j,» xl , C|J> 


1=0 


(a(a-j- (/-!)) 

l\ 


as x ^ 0. 


Applying Theorem 12.41 to and G^ yields their asymptotic expansions. 

Corollary 2.5. For iVyeR f/ie following asymptotic expansions hold as z ^ oo. 


yu) = z ‘e 2 E 


G"(z) = z'V y 


M-l 

E 

s=0 

M-l 

E 

s=0 


E %r(r + l/2) y bk,NCi, : 


j+r=s 


k+l-r 


y a,T(r + 3/2) y dk,NCl, : 


where aj are coefficients in the asymptotic expansion of K u w, 

' k -1 


K u (x) = x l ! 2 e x 


^ ajx i + 0(x k ) 
L i=o 


5 \ / ~ 


Fnt 1 (4^-(2/-i)2) 
2 8*j! 


= .9,v^(0)/j! are coefficients in the Taylor series about x = 0 of 


9n(x) = 


(x + l) w 
px + 2 


(2.14) 


(2.15) 


(2.16) 


+ 0(z~ 1 - M e~ z ), 

(2.17) 

+ 0(z~ 2 ~ M e~ z ), 

(2.18) 


(2.19) 


( 2 . 20 ) 
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dj,N = hx\o)/j\ are the coefficients in the Taylor series about x = 0 of 


and 


c i,j ~ c i,j- 1/2 


h N (x) = Vx + 2(x + 1) N , 

(-l/2-j)...(-l/2 -j-(l-l)) 


m\ 


are the coefficients in the Taylor series of (x + 1) '/ 2 J about x - 0. 


( 2 . 21 ) 

( 2 . 22 ) 


3 Recursion Relations 

Fjf(z) is analytic on 9t(z) > 0 and so it suffices to prove Theorem 11.21 for z e M + . z will be 
restricted to this domain from here on. The decay of the integrand, uniform in z, allows us 
to differentiate Fff (z) by differentiating under the integral. This will yield a pair of recursive 
formulas for Ff r (z). 

The derivations of the recursion relations do not depend on many of the details of the 
family of integrals in question, so we temporarily generalize to the family of integrals 


Hff{z) = f r(t) N C u (zr(t))w(t)dt 

Jti 


(3.1) 


where the only necessary assumptions are that the integrals are defined in the L 1 sense, 
differentiation commutes with the integral, and C(x) is a solution of the Bessel (upper sign) 
or modified Bessel (lower sign) equations 


~ 2 d 2 C u (z) dC u (z) _ . 2 2 , 


dz 2 

and satisfies the recursion relation 


d 


+ z- 


dz 


= (ir ^z 2 )C u (z) 


—[z- v C v {z)] = -z- v C v+l (z). 
dz 


(3.2) 


(3.3) 


These conditions are satisfied by Fff and G{f with C v = K v and taking the lower sign. 
Using Eq. (13.21) and differentiating under the integral gives a recursion relation for Hjf. 


Lemma 3.1. 


ttN+2/ .a 1 dH?(z) V* N ' 

H - + - * HAz) , 


(3.4) 


This reduces the task of computing Fff for non-negative integer N to the base cases 
IV = 0 and IV = 1, after which the recursion can be used to compute the result for higher N. 

Similarly, using Eq. (13.3ft and differentiating under the integral results in a second recur¬ 
sion relation. 


Lemma 3.2. 




(3.5) 
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Using Eq. (I1.5[) the base cases for Eq. (11.11) can be evaluated. 

Lemma 3.3. 

KO) = FjO) = lA' £ 5 .( 2 /2)A' 2 jl ( 2 /2). (3.6) 

Of particular interest are the cases IV = 0 and v odd or N = 1 and v even. In these cases 
the results involve modified Bessel functions of half odd integer order, which reduce to hnite 
linear combinations of the terms e~ z z~ l F~ n . In particular, we will need the following result. 

Lemma 3.4. K ±l / 2 (z) = (^)^ 2 z~ 1 ! 2 e~ z and for n- 1,2,..., 


Ah 


(n+ 1 / 2 ) (z) = * 1/2 6 Z |ri(2(ri-i) + l)(z n + z {n 1} ) + A n _ 2 (l/z)j (3.7) 


where R n ~ 2 is a polynomial of degree at most n - 2 . 
Proof. This is a consequence of 


FT _ z A {n + k)\z k 
±(n+ i/2)(") V 2 / 2 2 k k\(n - k)\ 


where n is a non-negative integer [T]. 

Combining Eq. (13.71) with Eq. (13. 6 j) implies 
Lemma 3.5. Ff(z) = § z~ l e~ z and for n > 0, 
2 2n 7r Y\™ll(2(n - j) + l ) 2 


K + i(z) =- 


z -l e -z ^-2n + z -(2n-l) + J R 2 „_ 2 (l/^)) , 


Ao(^) - VV*, A 2 1 (z) = ^rV^l + 2^), 


and for n > 1, 


Ait.') =- 77 ^ 7—77 - z l e (- ,2 " 11 + 2 lA * ' + fl 2 „- 3 (l/ 2 )) ■ 


(3,8) 

□ 


(3.9) 

(3.10) 

(3.11) 


2 ( 2 n-l) 

In the above formulas, A 2n _ 2 and A 2n _ 3 denote polynomials of degrees at most 2 n - 2 and 
2 n - 3 respectively. 

Using Lemma [3.51 and inducting via the recursion relation Eq. (13. 4 p yields the following. 

Lemma 3.6. For N odd and v even or vice versa, 

F?(z) = z- 1 e- z Q N!U (l/z) (3.12) 

where the Qn,v are polynomials whose degrees satisfy the following bounds for non-negative 
integers m and n, 


degQ 2 m, 2 n+i < 2 (n + m), degQ 2 m+i,o < 2 m, (3.13) 

deg Q2m+i,2n < 2 (n + m) - 1 for n > 0. (3.14) 
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Proof. Starting with the degrees of the base cases from Lemma 13.51 the recursion relation 
Eq. (13. 5 p implies that the degree increments by at most 2 every time N is increased by 2. □ 

In the following section we show that there are additional cancellations which cause the 
Qj t k s to have less than the maximum possible degree indicated here. 

Before moving on, note that another family of integrals of the form Eq. (13.11) for which 
similar results can be obtained using the recursion relation Eq. (13.41) is 

Hfi(z) - [ cos N (d)J„(zcos(6))d9, (3.15) 

Jo 

where the base cases are established via the following identity, found in Watson 5-43 (1) pLi]. 
Lemma 3.7. Let 93(yU + is) > -1. Then 


2 

71 


J fl+u (2zcos(6)) cos ((// - v)6)d6 = J ll {z)J v {z). 


This direction will not be pursued further here. 


(3.16) 


4 Proof of the Main Theorem. 

A proof of Theorem 11.21 will be given in this section. The main step involves using the 
recursion relations derived in the prior section to prove the following bound on polynomial 
degrees. 


Lemma 4.1. Recall the definition of the polynomials Qn,u for N odd and v even or vice 
versa, Eq. \3.FA) . Their degrees satisfy 


degQ^J "" 1 !/Ar< " 

’ [N-l if N > v. 

(4.1) 

For N odd and v even or vice versa. 


Gu ( z ) = z~ 2 e~ z P NtU (l/z) 

(4.2) 

where Pn,» are polynomials with 


deg P NtV < deg Q Nm - 1. 

(4.3) 

Proof. The case v - 0 follows from the bound Eq. (|3.13p. v = 

1 is the first example where 


nontrivial cancellation occurs. Using Eq. (13. 4 p one finds 


Ff(z) =^z 1 e z , Ff{z) = ^z 1 e Z (1 + 1 /z). (4.4) 

Note that the order of the polynomial increases by one, not two as would be naively expected 
from Eq. (I3.13p . The recursive formula Eq. (13.4p implies that the degree can increment by 
















at most 2 for each subsequent increase in N by 2. Explicit computation shows that it does 
increase by 2, at least initially. This proves the lemma for v - 1. 

We will assume v > 1 for the remainder of this proof. Using the recursion Eq. (|3.4[) gives 

F^ +2 (z) =z~ 1 e~ z [((1 - v 2 )z^ 2 + z~ x + 1)Q NiV (1/z) 

+ (3z~z + 2 z- 2 )Q' Nv (1/z) + z-*Cr^(llz)\ (4.5) 

=z~ 1 e~ z Q N+2 ,u- 


Suppose deg Qn,v < v - 1 for some N. Note that the bounds Eq. (13. 13jl imply that this 
holds initially i.e. 

deg QN min ,u <v-l (4.6) 

where N m[n - 0 for v odd and N min = 1 for u even. 

Denoting the coefficients of a polynomial Q by Qi we therefore have 


q u n%, =(i - v 2 )Qn:1 +3 (Q' N ,,r 2 + (Qn,uY~ 3 

= [(1 - V 2 ) + 3(z/ - 1) + (u - l)(i/ - 2)]Q£j (4.7) 

=0 


and 


Qn + 2 ,„ -(1 - v 2 )QnI + QnI + 3(Q(v,,r 3 + 2(Q(v,+ (Qtv^)^' 4 

=(2i/-l)(Q£j-Q5$J) (4.8) 

where polynomial coefficients of negative degree are defined to be zero. 

Therefore Qn+ 2 ,u has degree at most v and if Q'ffl = Q'ffl then the degree is in fact still 
bounded by v - 1. This allows us to conclude the following. 

Lemma 4.2. For v > 1 let L u denote the first positive integer that is even if v is even and 
odd if v is odd such that u + Qf^-i v - Then 

deg Q NtV <v -1 for N < L u , degQ N ,„ < u + N - (L u + 1) for N > L v . (4.9) 

Therefore, if L v - v for v > 2 then the proof of Eq. (14.1 ft will be complete. This will be 
shown by induction on v. 

From Eq. (13.101) . L 2 - 2 so the base case is proven. Suppose L u - u for some v > 2. Using 
Eq. (13.51) gives 


Qn+ i,i/+i(1/^) - [{v + 1)- 1 + l]Qiv,i/(l/ z ) + 21 2 Q'n,v(^/ z )- (4-10) 


Hence, for N < u 


Q 


v 

iV+l,i/+l 


= (1 + ")Qn,1 + Qn,u + (Qn,vY~ 2 

Q u n l 


(4.11) 
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and 


v-2 

V 


Qn+1,v+1 ~( U + ^)Qn^u + QnI + iy 2 )Qn, 

={2v - 1)Q%* + Q£j. 


(4-12) 


Therefore 


This proves Q 


V 

N+l.v+1 


= Q v ~ l 


Qn+1,u+1 Qn+1,v+1 ~ (Qn,u Qn,u) 


(4.13) 

)%+1 * Qu~u+1- Note that for V even ! 

1o„ +1 = Qo'J+i, as the recursion misses this case. This 
proves L v - v for all v > 2 by induction. This concludes the proof of Eq. (14.11) . 

Next we evaluate G^(z ) = F^ +2 (z) - F^(z). Using Eq. (11.71) yields 


Lemma 13.51 is needed to show that 0 


n+i, v +i for N + 1 < v and Q v v 
nv - n 

0,i/+l’ 


( Mn+2,u 

E 

s=0 

M N ,V 

- E 

s=0 


E a„r(r + 1/2) E ^k,N+2^m,r. 


k+m-r 


E o„r(r+l/2) E bk,N^m,'i 


k+m=r 


The 5 = 0 term of the expression in parentheses is 

a o c o,or(l/2)(&o,jv+2 - b 0jN ) = 0. 
Therefore, another z _1 term can be factored out to give 

G* ( z) = z~ 2 e~ z P NtU (l/z) 


(4.14) 

(4.15) 

(4.16) 

(4.17) 


where Pm v is a polynomial with degree at most M N+2 iI/ -l. This, together with the uniqueness 
of asymptotic expansions, proves Eq. (j4.2|) and Eq. (14.3|) and completes the proof this lemma. 

□ 


We conclude by showing that Theorem 11.21 is a simple consequence of the above lemma, 
combined with the asymptotic expansions from Corollary 12.51 

The uniqueness of asymptotic expansions in the asymptotic sequence e z z~ n (4j implies 
that the coefficients in the polynomials Qn,v and Pn,v must agree with those in the asymptotic 
expansions Eq. (12.171) and Eq. (12.181) respectively. The bounds Eq. (14.ip and Eq. (14.31) on the 
polynomial degrees therefore imply Theorem 11.21 
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A Neutrino Scattering 

In this appendix we outline the connection of the integrals Eq. (11.31) to neutrino scattering 
that motivated this paper. The following is based on results given in [6j. 

Consider 2-body to 2-body reactions between particles in kinetic equilibrium at a common 
temperature T and with distinct fugacities T$, i = 1,2,3,4. The reaction rate for the reaction 
1 + 2 -»• 3 + 4 in the Boltzmann limit is given by 


R =(T 1 T 2 - T 3 T 4 )/t, (A.l) 

i / r =2^1,0 A "i(^/ T ) ds ’ ( a - 2 ) 

where (|A4| 2 ) is the squared matrix element of the interaction, and r,r',t(y ) depend on the 
particulars of the reaction. See OT for definitions and further physics background. 

The matrix element for electron-positron annihilation into neutrinos, e~ + e + ++ v e + u e , is 
given by 


(\M\ 2 )=?>2G 2 F 


(1 + 2 y) 


,(s + t- ml) 2 


+ (2rj) 


:0 2 e-t) 2 


+ y(l + 2r))m 2 e s 


(A.3) 


where m e is the electron mass, y = sin(0 U) ) 2 , 9 W is the Weinberg angle, and Gp is Fermi’s 
constant. For this process 


r - \Js- 4m 2 , r - \fs , t(y) - ~(2rry - r 2 - (r') 2 ), 
r 1 IfiC 2 

y_ i (|Al| 2 )(s, t(y))dy = s [(16 y 2 + 8y - 1 )m 2 e + (8rj 2 + Ay + l)s]. 


(A.4) 


Therefore, the rate constant in the Boltzmann limit corresponding to this process is 


1 /r e - 


e' T <r+V e Ve 


2 10 7T 5 

g 2 f t 


oo ' 1 

Is ~fs {I i (l- M | 2 )( s > t (y))^) K i(PVs)ds 

n oo ._ _ 

J ^ sy/s - 4 m 2 [(16?r + 8y - 1 )m 2 e + (8 rj 2 + Ay + l)s] Ki(/3y/s)ds. 


(A.5) 


3(2 6 )tt 5 

Making the change of variables zcosh(f) = y/s/T, z = 2m e /T, the integral becomes 


1 /r e - 




GIT* j 


(2tt) £ 


r oo 

J /(cosh(f))A'i(^cosh(f)) sinh 2 (t)d£, (A.6) 


f(u) - Au 5 + Bu 3 , A = ~(8y 2 + Ay + 1), B - -(Ay 2 + 2rj - 1/4). 

3 3 

The integral Eq. (1A.6|) is of the form considered in Eq. (II.3ft with both N and v odd. Theorem 
11.21 does not apply but Theorem 12.41 can still be used to obtain the asymptotic expansions 
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as z -> oo of the relevant integrals 


(A.7) 

(A.8) 


G\(z) (l + \z-' + f- 

G»(z) =^' 2 e-» (l + jz-‘ + ^ 


_o 165 , 
z + ——z a + 


16 


2 1875 , 

Z 2 + ——z 6 + 
lb 


0(j- 4 )). 


Using Lemma [1.11 along with the recursion relation Eq. (13.4)1 also gives the exact formulas 
Gf(z) =^z- 1 K 0 (z/2) 2 + 2z- 2 K 0 (z/2)K 1 (z/2) + Q + ^z' 1 K^z/2) 2 , (A.9) 

G\{z) = Q + bz~ 2 ^ z _1 K 0 (z/2) 2 + (^ + 48z“ 2 J z~ 2 K Q {zj2)Ki{zl2) 


(: 


+ 1CU 2 + 96^ 


-4 


A ,(-/_>/- 


(A-10) 
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